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Abstract 

We have investigated the physics potential of very long baseline 
experiments designed to measure — > Ve oscillation probabilities. 
The principles of our design are to tune the beam spectrum to the 
resonance energy for the matter effect, and to have the spectrum cut 
off rapidly above this energy. The matter effect amplifies the signal, 
and the cut-off suppresses backgrounds which feed-down from higher 
energy. The signal-to-noise ratio is potentially better than for any 
other conventional beam experiment. 

We find that a beam from Fermilab aimed at the Super-K detector 
has excellent sensitivity to sin^(2^i3) and the sign of AM^. If the mass 
hierarchy is inverted, the beam can be run in antineutrino mode with 
a similar signal-to-noise ratio, and event rate 55% as high as for the 
neutrino mode. 

Combining the Fermilab beam with the JHF-Kamioka proposal 
adds very complementary information. We find good sensitivity to 
maximal CP violation for values of sin^(2^i3) ranging from 0.001 to 
0.05. 



" fritzd@fnal.gov 



1 



1 Introduction 



Measurements of the atmospheric neutrino flux with the Super-Kamiokande 
detector have shown that the deficit of muon neutrinos depends on the zenith 
angle as expected for neutrino oscillations This result combined with 
many other experimental constraints provides strong evidence for Vr 
oscillations, with mixing > 0.88% (90% C.L.) and AM^ in the range (1 — 
5) X 10-3 eV2 §. 

Efforts to study these oscillations with accelerator-produced neutrinos are 
in progress, using z/^ beams produced from pion decay. The K2K collabo- 
ration has detected neutrino interactions in the Super-K detector using a 
neutrino beam from KEK 250 km distant. They report a deficit of events 
compared to the expectation with no oscillations [Q. The MINOS experi- 
ment is under construction, and will use an Iron/scintillator detector 730 km 
distant from a neutrino beam produced at Fermilab f^. These experiments 
are expected to confirm that the disappearance has the LjE dependence 
expected for oscillations, and improve the precision on AM^. 

Measurements of solar neutrino fluxes also provide evidence for neutrino 
oscillations, usually explained as z/g — oscillations with Am^ much smaller 
than the AM^ for atmospheric oscillations. The Large Mixing Angle (LMA) 
solution is strongly favored by the data, with Am^ in the range (3 — 25) x 10"^ 
eV^ HI- The KamLAND and Borexino experiments, currently under 
construction, are expected to test the LMA solution and better constrain the 
mixing parameters. 

The atmospheric and solar neutrino oscillations can be simultaneously 
explained with mixing in a three-flavor model: The the three known flavors 
of neutrinos, v^,, u^, and u^-, are mixtures of mass eigenstates i^i, and 
z/3. The data require a mass hierarchy, with |m3 — 7712! ^ \m2 — The 
hierarchy can be "normal", with ms > m2, or "inverted", with m2 > rris. We 
identify AM^ = m| — and Am^ = — ml. 

While z/g oscillations are mainly associated with solar oscillations, 

they can also occur at the atmospheric frequency if there is a z/3 component 
in the z/g eigenstate. Reactor experiments have placed limits on the size of 
this mixing, resulting in the limit sin^(26'i3) < 0.1 P, 0. If sin^(2^i3) > can 
be detected, the phenomenology of z/^ — z/g oscillations becomes very rich: 

• Interference between the solar and atmospheric oscillation amplitudes 
can result in large CP violation [l^ . 
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Figure 1: CP asymmetry at the peak of the atmospheric oscillation for the 
parameters of Table |l|. Solid line is for Am^ = 5.0 x 10~^ eV^, and dashed 
line is for Am^ = 2.0 x 10"^ eV^. 



Matter effects can modify the oscillation probability, depending on the 
sign of AM2 [O, P, p 



To study these effects on the oscillation probability, we use the paramet- 
rization of Ref. [jl4|, with the oscillation parameters of Table ^ as defaults. For 
example. Fig. |T] shows the possible CP asymmetry as a function of sin^(2^i3), 
at the peak of the atmospheric oscillations. While the mixing in solar os- 
cillations is greater than in atmospheric oscillations, the solar oscillations 
develop much more slowly, so the two amplitudes for — > z/g can be com- 
parable at atmospheric baselines. We note that maximal CP asymmetry is 
a possibility: There could be a large signal for neutrinos and no signal at all 
for antineutrinos, or vice-versa. 

Given the potential for profound insights into fundamental physics, most 
studies of future neutrino oscillation experiments with conventional beams 
have focused on z/„ u^, oscillations. A wide variety of neutrino energies and 



detector technologies have been considered in Refs. [15, 16, 17, 18 
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Table 1: Default oscillation parameters 



The potential of low energy neutrino beams combined with large water- 



cerenkov detectors has been explored for beams based at JHF [T^ , CERN pO 



and Fermilab ITH] and baselines of 100 to 300 km. These low energy ap- 
proaches have the advantage of very low backgrounds, and have excellent 
sensitivity to very small z/g appearance signals. They have the disadvantage 
that matter effects are too small to disentangle from CP violating effects, 
but are still large enough that ambiguities from the unknown sign of AM^ 
affect the ability of the experiment to establish CP violation pT |. 

The potential of higher energy neutrino beams with very long baseline 
(L > 7000 km) has been explored in Refs. |T3|, [T3]. The backgrounds in 



water cerenkov detectors are larger for higher neutrino energies. However, the 
i^ti signal at these baselines can be highly amplified by matter effects, 

improving the signal-to-background ratio [jll] , |12|, |13[. Also, the amplified 
signal depends mostly on atmospheric oscillations, with solar oscillations 
playing a much smaller role than with shorter baselines. 

From this summary of short- and long-baseline neutrino oscillation ex- 
periments, we see that they can play very complimentary roles. The short- 
baseline experiment measures a CP-violating combination of atmospheric 
and solar oscillations. The long-baseline experiment determines the sign of 
AM^, constrains the matter effects, and constrains the size of the atmo- 
spheric z/^ — > Vf, oscillation alone. 

We propose herein an optimized long-baseline experiment to complement 
a short-baseline experiment. For concreteness, we assume a target detector at 
the Kamioka site in Japan, where the Super-K detector already exists. The 
JHF-Kamioka proposal [jl9l uses the Super-K detector for the first phase, and 
for the second phase assumes construction of a new detector with 40 times 
more fiducial mass. We assume that this detector will be used as the target 
for two neutrino beams: One from JHF, with a baseline of 295 km, and one 
from Fermilab, with a baseline of 9300 km. The concept of our proposal 
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does not depend on the details of this choice, and could be adapted for other 
locations. The main principles of our proposal are: 

1 . Combine information from very short and very long baselines in order to 
determine the mass hierarchy and CP-violating oscillation parameters. 

2. Match the spectrum of the long-baseline neutrino beam to the energy at 
which matter effects produce the maximum amplification of the signal. 

3. Design the long-baseline neutrino beam to have an energy spectrum 
with a rapid cut-off above the signal region. Since most backgrounds 
feed-down from the neutrino energy to a lower visible energy, this re- 
duces the background in the signal region. 



2 The z/^ z/g probability at 9300 km 

We start in the "leading approximation" which parametrizes oscillations 
driven by the atmospheric AM^ and neglects oscillations driven by the solar 
Am^. The t'^ z/g probability can then be written as 

Piu, v^) = Rl sin^ e,, sin^ 2^13 sin^ A-267 AM^ M _ 

\ -Km rlju J 

The effects of interaction with matter are included in the parameter Rm- 
In vacuum, Rm = 1. For positive AM^, there is a resonant energy at 
which there is full mixing, with R^ = 1/ (sin ^23 sin 2^13). For antineutrinos, 
this resonance occurs for negative AM^. Since we're considering the case 
sin^2^i3 < 0.1, the oscillation length becomes very long, and for terrestrial 
baselines we can expand the probability to first order in L: 

P(z/^ ^z/e) = sin2 ^23 sin^ 2^13(^1. 267 AM^A^ . (2) 

Since the flux is proportional to l/L^, we find that the number of expected 
signal events is independent of distance. Since background rates decline with 
L^, longer baselines will provide better signal to noise. For negative AM^, 
the signal would be suppressed and unobservable. 

Fig. ^ shows -P(t';i i^e) as a function of neutrino energy for a baseline 
of 9300 km, AM^ = 3.5 x 10"^ eV^ and sin^ 2^13 = 0.01. We see that the 
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Figure 2: Oscillation probability 
as a function of energy, for L = 
9300 km, sin2(2^i3) = 0.01, and 
neglecting solar oscillations. Solid 
histogram is from a numerical cal- 
culation p2|, dashed line is from 
the analytical approximation of 
Ref. [|I^], with the matter parame- 
ter A set to 1.75 X 10"^ eVVGeV. 



Figure 3: Oscillation probability 
function of energy, for three 
different values of the atmospheric 
AM^ and parameters of Fig. |[ 



resonant energy is in the range 8 to 10 GeV. At this energy, the oscillation 
probability is 20 times the maximum probability occurring in vacuum. 

Fig. also compares a numerical calculation of the exact theory in- 
corporating the density variation of the Earth, with an analytical approxima- 
tion 1 14] using an average density. The agreement is excellent in the energy 
range of interest, and we use the analytical approximation for the rest of this 
paper. Refs. |Tl], [1^ have also shown that using an average density is a good 
approximation for similar parameters. 

Fig. ^ shows P{iy^ — > t'e) for different values of AM^. The resonant 
energy is proportional to AM^, while the maximum oscillation probability is 
constant. 

The analytical approximation of Ref. jl^ can also be used to include 
non-leading terms. The amplification occurs mainly for the leading terms, so 
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the oscillation probability for a 9300 km baseline is dominated by the leading 
terms much more than for shorter baselines. 



3 The neutrino beam 

Our starting point for the neutrino beamline design is the NuMI beamline at 
Fermilab, which is currently under construction. In the NuMI beamline, 120 
GeV protons from the Main Injector collide with a target which is designed 
for an average power of 0.4 MW. Pions produced in the target are focused 
and charge-selected in a double horn system, and then decay in a 675 m 
decay tunnel. This produces a beam consisting of mainly z/^. Different 
configurations of the target and horns produce spectra peaking at different 
energies. There is a small z/g contamination arising from the decay of kaons 
produced in the target, and the decay of muons produced in pion decays. 
The level of this contamination varies from 0.6 to 1%, depending on the 
configuration. 

The accelerator used to inject protons into the Main Injector is the 
Booster ring, which accelerates protons up to 8 GeV. Proton drivers with 
power of 1 MW or more are currently being designed p3|, . With a proton 
driver replacing the Booster, the Main Injector would be capable of deliver- 
ing 2 MW of proton power to target. In the case of the linac option for the 
proton driver, it may be possible to deliver this power at any proton energy 



less than 120 GeV \24 



To aim a neutrino beam from Fermilab to the Kamioka site, the decay 
tunnel must be angled downwards 46° from the horizontal. At the Fermilab 
site, the decay tunnel must fit within a 200 m vertical depth in order to be 
within rock favorable for construction and to avoid a water aquifer. This 
limits the length of the decay tunnel, and reduces the efficiency for pions to 
decay and produce neutrinos. The length is further reduced by the vertical 
space needed to bend the protons down to the target, and the space needed 
for the target and focusing, and beam dump at the end of the tunnel. 

As discussed in Sect. ^ most backgrounds for our measurement feed-down 
from the neutrino energy to a lower visible energy, while for the signal the 
entire neutrino energy is visible. Backgrounds in the signal region tend to 
arise from neutrinos with higher energy. Therefore, our goal for the neutrino 
beamline is to produce a spectrum which peaks at the energy for the maxi- 
mum oscillation probability in Fig. 13, and falls off quickly at higher energy. 
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Figure 4: Solid line: Spectrum of interacting CC events at the Kamioka 
site, assuming no oscillations. Dashed line includes effects of z/^ disappear- 
ance from oscillations. We note that the second disappearance peak happens 
to occur at the peak of the spectrum. 



A saw-tooth shape would be ideal. The lack of a high-energy tail might also 
simplify the radiation shielding requirements relative to NuMI. 

The beamline will need to be tuned for the actual value of AM^ when it 
is better known; for now we assume AM^ = 3.5 x 10~^ eV^, and aim for a 
peak energy of 9 GeV. 

We have used the NuMI fast beamline simulation to produce a candidate 
beamline design. For the following, we apply a scaling factor of 1.25 needed 
to bring the fast simulation into agreement with the full simulation of the 



medium- and high-energy NuMI configurations [25|. Starting with a standard 



NuMI configuration, we have made the following changes: 

1. Use 60 GeV protons on target instead of 120 GeV. This reduces the cost 
of the bending magnets, and allows the bend to occur in less vertical 
space. We assume a field of 6 T for the bending magnets. 
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2. We assume the accelerator complex is upgraded to provide 2 MW of 
proton power from the Main Injector. 



3. Change the target and horn configuration to produce a spectrum peak- 
ing at 10 GeV. 

4. Add a small dipole bend after the horns, to bend 20 GeV pions by 0.5°. 
This allows the higher energy pions to be filtered out, and reduces some 
of the z/e component by filtering out neutral kaons. 

The resulting spectrum is shown in Fig. ^, and has the basic properties we 
were aiming for. The spectrum including oscillation effects is also shown. 
In an interesting coincidence, the 2nd disappearance peak happens to occur 
at the peak of the spectrum. The length of the pion decay tunnel is 210 m, 
and the decay efficiency is 37% what it would be for a NuMI length tunnel. 

The event rates do not depend critically on the specifics of the design. For 
example, if we need to use 120 GeV protons to get the full beam power, the 
pion decay efficiency relative to NuMI decreases only to 35%. Or, if we use 
2 T bending magnets and 60 GeV protons, the relative efficiency decreases 
to 33%. 



Ref. has also presented beamline simulations producing energy spec- 
tra with sawtooth shapes, using quadrupole magnets instead of horns to focus 
the pions. 



3.1 Flux normalization 

Since the use of space for a near detector would further reduce the efficiency 
of the beamline, we have assumed that there will be no near detector and the 
flux can be normalized using data in the far detector. As shown in Fig. ^, we 
expect a large rate of z/^ charged-current interactions in the far detector, and 
this can in principle be used to normalize the flux. The systematic limitations 
of this procedure will have to be studied with a detailed simulation. It's 
possible that neutral-current interactions or 1/^- charged-current interactions 
might also be used to help measure the flux. 



3.2 Running with antineutrino beams 

The neutrino beamline can be converted into an antineutrino beamline by 
reversing the polarity of the focusing horns and bending magnet. Since the 
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antineutrino cross-section is smaller than the neutrino cross-section, and vr" 
production in the proton target is suppressed relative to 7r+ production, event 
rates in the far detector are typically 1/3 of those in neutrino beams. The 
intrinsic z/g and Ve backgrounds tend to be worse as well. These problems are 
mitigated by the following effects: 

• As discussed in Sect. H, the efficiency of reconstructing z/g events is 
higher than for z/g events. After reconstruction and event selection 
cuts, the antineutrino rate will be about 55% of the neutrino rate. 

• The dipole bend helps remove neutral kaons and in the pion decay 
tunnel, removing the v^. component of the beam. The remaining 
component will be comparable to the v^. component of the neutrino 
beam. 

4 The z/g signal 

The Vf, signal is detected from charged- current (CC) interactions producing 
an electron along with a hadronic shower with energy -Ehad- The visible 
energy of the event, i?vis, (the energy not carried away by neutrinos) is equal 
to the neutrino energy. The differential cross-section in terms of the variable 
y = EYiat^l Ey, normalized to 1, can be approximately parametrized by pPj ]: 



We will model the efficiency of detecting a Ue CC event by assuming 
the fiducial mass of the detector for the target mass, and by requiring the 
production of an electron above some threshold E. The probability of a Ve 
with energy E^ generating an electron above a threshold E is: 




(3) 



For antineutrinos, the parametrization is: 
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Figure 5: Solid line: Spectrum 
of interacting v^, CC events at 
the Kamioka site, assuming the 
oscillation parameters of Fig. H. 
Dashed line includes effects of cuts 
described in Sect. 0. 



Figure 6: Solid line: Spectrum 
of interacting v^. CC events at 
the Kamioka site, assuming the 
oscillation parameters of Fig. ^. 
Dashed line includes effects of cuts 
described in Sect. ^. The rate after 
cuts is 55% of the rate for v^, CC 
events. 



For antineutrinos, the parametrization is: 



Vc 



15 / 8 
y I 15 



1 E 
5^ 




(6) 



For neutrinos, we choose the following cuts for the threshold E: 

1. E> 4.5 GeV. 

2. E/^vis > 0.45. 

Fig. ^ shows the z/g CC energy spectrum at the Kamioka site, assuming the 
neutrino beam described in Sect. ^, and the neutrino oscillation parameters 
of Fig. 1^. The dashed line of Fig. ^ shows the effect of the above cuts. The 
average efficiency of these cuts is 50%. With these cuts, we find 0.21 events 
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in the signal region, per kton-year of exposure, where the signal region is 
defined as energy between 7 and 10.5 GeV. 

The ^/-distribution for antineutrinos is more peaked towards zero than 
for neutrinos, as described by Eqs. and (^. Therefore, for the same NC 
background fraction, z/g CC events can be selected with more efficient cuts. 
For neutrinos, we choose the following cuts for the threshold E: 

1. E> 3.3 GeV. 

2. E/E^is > 0.33. 

Fig. P shows the Ue CC energy spectrum at the Kamioka site, assuming the 
neutrino beam described in Sect. ^, and the neutrino oscillation parameters 
of Fig. ^ with negative AM^. The dashed line of Fig. ^ shows the effect of 
the above cuts. The average efficiency of these cuts is 80%, producing an 
antineutrino signal rate 55% of the neutrino rate. 



5 Backgrounds 

5.1 Backgrounds from neutral-current events 

Neutral pions are produced in the hadronic showers present in all neutrino in- 
teractions with nuclei. At high energies, the photons from the tt^ '-fj decay 
produce a single electromagnetic shower, which in water cerenkov detectors 
is indistinguishable from that of an electron. Neutral-current (NC) events 
with an energetic 7r° constitute a major background to the z/g appearance 
signal. The visible energy of these events is given by yE,y, so they feed-down 
from the neutrino energy to a lower visible energy. 

We have studied the production of neutral pions with the LEPTO pro- 



gram which generates deep inelastic neutrino interactions and models the 
fragmentation of the hadronic component. There are several efforts underway 
to model resonance and quasi-elastic production, which we have not included, 
and which might account for roughly 15% of the cross-section in the energy 



range we are interested in ||2^. However, we do expect the backgrounds to 
come more from high-?/ events, or from feed-down from higher-energy neu- 
trinos, where the fragmentation assumption is more applicable. 

We find that Eq. (|^) and Eq. (^) agree well with LEPTO-generated data 
for charged-current and neutral-current events. 
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Figure 7: Integral plots of the E{7r^)/E{h.a.d) distribution in neutrino inter- 
actions. The different histograms are for different ranges of hadronic energy: 
Solid: 20 to 25 GeV. Dashed: 15 to 20 GeV. Dotted: 10 to 15 GeV. Dot- 
Dashed: 5 to 10 GeV. The curve is our parametrization of this quantity. 



The 7r° fragmentation functions are roughly independent of hadronic 
energy. This is illustrated in Fig. which shows integral plots of the 
E{7r^) I i?(had) distribution. The probability for a hadronic shower to produce 
a with energy greater than a fraction x of i?(had) can be parametrized 
by: 

= (0.49) - (0.96)x + (0.47)xl (7) 

For a given neutrino energy, we can use Eq. |^ to obtain the visible energy 
spectrum in NC interactions. For a given visible energy, we can then use 
Eq. 1^ to calculate the probability that there will be a 7r° above some cut-off 
energy. Thus, we can convolute an input neutrino spectrum into an output 
NC background spectrum. 

Fig. I shows the NC background to the z/g signal, estimated with this 
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Figure 8: Solid line: Spectrum of 
z/e events with cuts from Fig. ^. 
Dashed line: Distribution of back- 
ground from NC events. 



Figure 9: Solid line: Spectrum of 
Ve events with cuts from Fig. ^. 
Dashed line: Distribution of back- 
ground from NC events. 



procedure with the cuts of Sect. We find a background of 0.047 events 
per kton-year in the signal region, corresponding to fs = 2.1%, where fs is 
the background fraction relative to the expected signal after cuts for a beam 
with 100% z/g [|r3]. Fig. ^ shows the same distributions for the z/e signal. We 
obtain the same background fraction with a signal rate 55% of the Vf. rate. 

Since the direction of the incoming neutrino is known, and the outgoing 
neutrino can carry away transverse energy, the direction of the observed 
energy flow might also be used to further reduce NC background. A detailed 
simulation will be needed to investigate this possibility. 



5.2 Background from charged- current interactions 

In this analysis, we assume that all z/^ CC events can be rejected due to the 
presence of a muon, and only consider the vr" background from NC events. 
For high-?/ events however, the muon will have low energy and can easily be 
missed, while the hadronic shower has most of the energy of the neutrino. 
Neutral pions from these hadronic showers will generate backgrounds the 
same as for neutral currents. Fortunately, as illustrated in Fig. ^, the z/^ flux 
is very low in the signal region. This is true for other values of AM^ as well. 
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So, we expect this source of background to be negligible compared to the 
background from neutral current events. 



5.3 Background from Vt charged-current interactions 

At the peak energy for the signal, most of the t'^ have oscillated into Vri so 
Vr CC interactions are a potentially serious source of background. However, 
at this energy the Pr CC cross-section is only ^ 8% of the Pe CC cross- 
section. The decay mode r — > euu produces a real electron background, but 
the branching ratio is only 18%. Neutral pions in Vj. CC events can also 
produce backgrounds. Since all these events have an escaping neutrino from 
the r decay, these backgrounds feed-down to lower visible energy. As long 
as we keep the high energy tail of the neutrino beam suppressed, we expect 
this source of background to be small. However, the cuts needed to suppress 
this background need more study with a detailed simulation. 



5.4 Intrinsic z/g backgrounds 

As discussed in Sect. ^ the neutrino beam has a small amount of Ve contam- 
ination. The dipole bend in the beamline can help reduce this background 
by eliminating neutral kaons in the beamline. Also, the Ue energy spectrum 
is very broad, while the signal region is narrow. The intrinsic Ve background 
as a function of the energy width of the beam has been studied in Ref. ]TB 



showing that it is greatly reduced for narrow beams. Therefore, we expect 
this background also to be very small. 



6 Comparison with the JHF-Kamioka project 

The first stage of the JHF-Kamoika project uses a 0.75 MW proton driver 
to generate a neutrino beam aimed at the Super-Kamiokande detector. In 
phase 2, the proton power is increased to 4 MW, and a new water cerenkov 
detector, called Hyper-K, is built with 40 times more fiducial mass. 

With very similar mixing parameters as those used in Sect. ^, the JHF- 
Kamioka proposal obtains 0.11 events per kton-year in phase 1, compared to 
0.21 for our proposal. The event rate for the Fermilab beam is 200% of the 
rate for the phase 1 JHF beam, and 40% of the rate for the phase 2 JHF 
beam. 
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The JHF-Kamioka project achieves a very low background rate: = 
0.51%, compared to J'b = 2.5% for the Fermilab beam, where J'b is the 
background fraction relative to the expected signal after cuts for a beam 



with 100% Ve [|15l. However, the signal from the Fermilab beam is amplified 
by a factor of 20, so the effective background for comparison with the JHF 
beam is only 0.13%. 

The signal region for the Fermilab beam is fairly narrow, between 7 and 
10.5 GeV. The signal region for the JHF beam is defined as energy between 
0.4 and 1.2 GeV. The average oscillation probability over the JHF spectrum 
is only 56% of the peak probability. This further accentuates the background 
advantage of the Fermilab beam. 



7 Combined analysis of data from Fermilab 
and JHF beams 

In the previous section, we showed that the Fermilab and JHF beams would 
have comparable performance. The JHF phase 2 beam has an advantage in 
rate, while the Fermilab beam has an advantage in background. However, 
the real advantage comes in combining information from the two beams. The 
first step will be the detection of u^. oscillations and determination of 

the sign of AM^. The second step will be constraints on sub-leading effects 
such as CP violation. 



7.1 Detecting oscillations and the sign of AM^ 

Fig. ^ shows our estimates of the reach in sin^(2^i3) for JHF and Fermilab 
beams with the Super-K detector, assuming positive AM^, and atmospheric 
oscillations only. The JHF curves are in reasonable agreement with curves 
in the JHF-Kamioka proposal [0. If a signal appears with the JHF beam, 
it should be confirmed by the Fermilab beam within one year of running, 
otherwise AM^ must be negative. It should then be possible to find the 
signal with the Fermilab beam in antineutrino mode. The performance of 
the Fermilab beam is excellent: 5 years of running yields a 3a reach of 0.003 
in sin2(2^i3). 

Fig. |ll] shows a similar plot for the Hyper-K detector. A 5 year run 
should produce a 3a reach of 0.001 or better, with a strong dependence on 
the background level and systematic uncertainty. 
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Figure 10: Reach in sin^(20i3) 
with Super-K detector. Upper 
curves are for JHF phase 1 beam, 
lower curves for Fermilab beam. 
Sohd curves assume 10% sys- 
tematic uncertainty on the back- 
ground, dashed curves assume no 
systematic uncertainty. We as- 
sume positive AM^, and atmo- 
spheric oscillations only. 



Figure 11: Reach in sin^(2^i3) 
with Hyper-K detector. Upper 
curves are for JHF phase 2 beam, 
lower curves for Fermilab beam. 
Solid curves assume 10% sys- 
tematic uncertainty on the back- 
ground, dashed curves assume no 
systematic uncertainty. We as- 
sume positive AM^, and atmo- 
spheric oscillations only. 



It is interesting that if sin^(2^i3) is small enough for Eq. (0) to be valid, 
the long-baseline oscillation probability does not depend explicitly on the 
size of the matter effect for neutrinos at the resonant energy. If the peak 
energy for the matter resonance can be determined, the theory of matter 
effects can be tested and systematic uncertainties from matter effects can be 



kept small 13 



7.2 Constraints on the CP violating phase 

The JHF-Kamioka phase 2 program assumes 2 years of neutrino running, and 
6 years of antineutrino running. They have assumed the background will be 
very well constrained, and for our analysis we have assumed a 2% systematic 
uncertainty on the background. Except for the highest allowed values of 
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ye- 




sin^(2^i3) 




Beam 


ars 


0.01 


0.05 


0.001 


Fermi 


2 


0.104 ±0.0060 


0.501 ±0.012 


0.012 ±0.0034 


Fermi 


8 


0.104 ±0.0032 


0.501 ±0.0061 


0.012 ±0.0021 


JHF V 


2 


(7.9 ±0.37) ■ 10-3 


(3.2 ±0.052) ■ 10-2 


(1.5 ±0.32) ■ 10-3 


JHF V 


8 


(7.9 ±0.25) ■ 10-3 


(3.2 ±0.031) ■ 10-2 


(1.5 ±0.24) ■ 10-3 


JHF V 


6 


(2.8 ±0.33) ■ 10-3 


(1.8 ±0.044) ■ 10-2 


(0.0 ±0.31) ■ 10-3 



Table 2: Expected oscillation probability measurements for the various beams 
and input values. 



sin2(26'i3), it will be necessary to either find a way to further reduce the 
background, or to determine the background to this level. For the Fermilab 
beam, we have assumed a 5% systematic uncertainty on the background, and 
an 8 year run. 

We have neglected the systematic uncertainty on the fiux prediction, but 
note that for the highest allowed values of sin^ (26*13), it will need to be 
constrained to 1%. 

For the JHF beam, we have assumed the signal to background ratio in 
the antineutrino beam is the same as for the neutrino beam, and the 6 year 
exposure for the antineutrino beam is equivalent to 2 years with the neutrino 
beam. Without a dipole bend, it is likely that the backgrounds with the 
antineutrino beam will in fact be much worse. As an alternative, we have 
considered running for 8 years with the JHF neutrino beam. 

We have assumed the parameter values in Table |I], and sin2(26'i3) = 0.01, 
0.05, or 0.001. We have not included the effect of uncertainties in parameters 
other than sin2(2^i3) and b. 

We have assumed that the JHF beam measures oscillation probabilities 
at 0.8 GeV, and the Fermilab beam at 9.0 GeV, although in reality a more 
sophisticated energy-dependent analysis will be performed. The expected 
oscillation probability measurements are shown in Table ^. 

Constraints in the ^\Y?{2Q\-i)-b plane from various combinations of these 



measurements are shown in Figs. through ^ For each input value of 



sin2(2^i3), maximal CP violation can be detected at 3(T or better, demon- 
strating a sensitivity over an impressive range of parameter space. This is 
possible since, as illustrated in Fig. |T], the maximum CP asymmetry is larger 
if sin2(2^i3) is smaller, so the larger asymmetry provides statistical com- 
pensation for the smaller signal. Over most of this range, 2 years with the 
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Fermilab and JHF beams combined is enough to reach the 3a criterion. The 
constraints are improved with an 8 year run during which the JHF beam 
runs in antineutrino mode for 6 years. If the JHF antineutrino beam does 
not provide adequate performance, an 8 year run with Fermilab and JHF 
neutrino beams provides constraints almost as effective. 

8 Summary and conclusions 

We have investigated the physics potential of very long baseline experiments 
designed to measure oscillation probabilities. The principles of our 

design are to tune the beam spectrum to the resonance energy for the matter 
effect, and to have the spectrum cut off rapidly above this energy. The 
matter effect amplifies the signal, and the cut-off suppresses backgrounds 
which feed-down from higher energy. The signal-to-noise ratio is potentially 
better than for any other conventional beam experiment. 

The measurements from long- and short-baseline experiments are very 
complementary. The short-baseline experiment measures a CP-violating 
combination of atmospheric and solar oscillations. The relatively small mat- 
ter effects can still confuse the interpretation of the measurements. At long 
baselines, solar oscillations play a relatively much smaller role. The long- 
baseline experiment determines the sign of AM^, constrains the matter ef- 
fects, and constrains the size of the atmospheric z/^ oscillation alone. 

In particular, we have investigated the capabilities of a neutrino beam 
from Fermilab aimed at the Super-K detector in Japan, 9300 km distant. At 
this baseline, neutrinos at the resonance energy have an oscillation proba- 
bility 20 times the maximum in vacuum. We have used a simple model to 
estimate the signal efficiency and background, and it will be important to 
check this with a full simulation. For a normal mass hierarchy, a five year 
run can detect sin^(2^i3) at 3a for values as low as 0.003. If the mass hier- 
archy is inverted, the beam can be run in antineutrino mode, with a similar 
signal-to-noise ratio, and event rate 55% as high as for the neutrino mode. 

In an interesting coincidence for the 9300 km baseline, the second peak 
of maximum z/^ disappearance occurs at the same energy as the peak for 
appearance. Thus, interesting measurements of the parameters of disap- 
pearance may also be possible. 

Phase 2 of the JHF-Kamioka proposal assumes the construction of the 
Hyper-K detector with 40 times the fiducial mass of the Super-K detector. 
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Having this detector serve as the target for beams from both JHF and Fermi- 
lab is a powerful way of placing constraints in the sin^ {29i-^)-6 plane. We have 
investigated values of sin^(26'i3) ranging from 0.001 to 0.05, and find that an 
8 year run can establish maximal CP violation at 3a or better throughout 
this range. If the JHF antineutrino beam has poor performance, using Fer- 
milab and JHF neutrino beams alone is almost as effective in constraining 
S. 
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Input values: sin^(26'i3)=0.01, 6 = 7i/2. Contours are 1, 2, 3a 




0.006 0.008 0.01 0.012 0.014 
sin (20^3) 



0.006 0.008 0.01 0.012 0.014 
sin (20^3) 



Figure 12: After 2 years with Fer- 
milab beam. 



Figure 13: After 8 years with Fer- 
milab beam. 
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Input values: sin^(26'i3)=0.01, 6 = 7t/2. Contours are 1, 2, 3a 



3 - 




■■a - 



to 

-1 - 



-2 - 



■■a 
to 

-1 



0.006 0.008 0.01 0.012 0.014 
sin (20^3) 

Figure 16: After 2 years with Fer- 
milab and JHF beams. 



3 F 




0.006 0.008 0.01 0.012 0.014 
sin'(26j3) 

Figure 18: After 2 years with JHF 
beam and 6 years with JHF an- 
tineutrino beam. 



0.006 0.008 0.01 0.012 0.014 
sin (20^3) 

Figure 17: After 8 years with Fer- 
milab and JHF beams. 



3 F 




-2 - 



0.006 0.008 0.01 0.012 0.014 
sin'(26j3) 

Figure 19: After 8 years with Fer- 
milab, 2 years with JHF, and 6 
years with antineutrino beams. 
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Input values: sin^(26'i3)=0.05, 6 = 7i/2. Contours are 1, 2, 3cr 




0.04 0.045 0.05 0.055 0.06 0.04 0.045 0.05 

sin^(2ej3) sm\2Q^^) 



0.055 



0.06 



Figure 20: After 2 years with Fer- Figure 21: After 8 years with Fer- 
milab beam. milab beam. 




0.04 0.045 



0.05 



sin (2ej3) 



0.055 0.06 0.04 0.045 0.05 



0.055 



0.06 



Figure 22: After 2 years with JHF Figure 23: After 6 years with JHF 
beam. antineutrino beam. 
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Input values: sm^(26'i3)=0.05, 6 = 7i/2. Contours are 1, 2, 3a 



3 - 




-2 - 

^ ^ ^ 

0.04 0.045 0.05 0.055 0.06 
sin^(2ej3) 

Figure 24: After 2 years with Fer- 
milab and JHF beams. 



3 - 




-1 - 



-2 - 



0.04 0.045 0.05 0.055 0.06 
siii^(2ej3) 

Figure 26: After 2 years with JHF 
beam and 6 years with JHF an- 
tineutrino beam. 



3 - 




-2 - 

0.04 0.045 0.05 0.055 0.06 
sin^(2ej3) 

Figure 25: After 8 years with Fer- 
milab and JHF beams. 



3 - 




-2 - 



0.04 0.045 0.05 0.055 0.06 
sin^(2ej3) 

Figure 27: After 8 years with Fer- 
milab, 2 years with JHF, and 6 
years with antineutrino beams. 
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Input values: sin^(26'i3)=0.001, 6 = 7i/2. Contours are 1, 2, 3a 




0.001 0.002 0.003 0.004 
sin\20.,) 



0.001 0.002 0.003 0.004 
sin^(20.,) 



Figure 28: After 2 years with Fer- Figure 29: After 8 years with Fer- 
milab beam. milab beam. 




0.001 0.002 0.003 0.004 
sin^(20.^) 



0.001 0.002 0.003 0.004 
sin^(20.,) 



Figure 30: After 2 years with JHF Figure 31: After 6 years with JHF 



beam. 



antineutrino beam. 
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Input values: sin^(26'i3)=0.001, 6 = 7i/2. Contours are 1, 2, 3a 




sin^(20j3) 



0.002 0.003 0.004 
sin^(20.,) 



Figure 32: After 2 years with Fer- Figure 33: After 8 years with Fer- 
milab and JHF beams. milab and JHF beams. 




0.001 0.002 0.003 0.004 

sin^(20j3) 

Figure 34: After 2 years with JHF 
beam and 6 years with JHF an- 
tineutrino beam. 



0.001 0.002 0.003 0.004 

sin^(20j3) 

Figure 35: After 8 years with Fer- 
milab, 2 years with JHF, and 6 
years with antineutrino beams. 
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